Cross-Kerr interaction in a four-level atomic system 
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We derive the form of the cross-Kerr interaction in a four-level atomic system in the N- 
configuration. We use time-independent perturbation theory to calculate the eigenenergies and 
eigenstates of the Schrodinger equation for the system. The system is considered as a perturbation 
of a Raman resonant three-level lambda scheme for which exact solutions are known. We show 
that within the strong control field limit the cross-Kerr interaction can arise between two weak 
probe fields. The strength of this nonlinear coupling is several orders of magnitude larger than that 
achievable using optical fibres. 
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I. INTRODUCTION 

The cross-Kerr interaction involves controlling the re- 
fractive index experienced by one electromagnetic field by 
the intensity of another. This nonlinear coupling between 
the field modes is vital in several quantum information 
protocols [1, 0, @1 has found applications in nonde- 
molition photon measurement [l| and entanglement con- 
centration [1, 

However, achieving a strong nonlinear interaction 
whilst avoiding self-phase modulation or large absorp- 
tion has proved challenging. Conventional methods in- 
volve either four-wave mixing in three-level atomic sys- 
tems [1] or rely on the weak third-order susceptibility 
experienced in microstructured optical fibres H. It has 
been suggested by Schmidt and Imamoglu [1(J| that a 
much stronger interaction is generated in the four-level 
atomic system (Fig. [T]). The experimental feasibility of 
achieving a large cross-Kerr interaction for continuous- 
wave fields has been demonstrated However it is 
only recently that a suitable method was suggested for 
group- velocity matched pulses [13 , as might be necessary 
for the construction of quantum logic gates. 

The Hamiltonian for a cross-Kerr interaction [?[ be- 
tween two fields J7a and fie is given by 



H — hKflnflr 



(1) 



where K is the coupling strength and fia and he are 
photon-number operators. When acting on a state con- 
sisting of two modes, both of which are in number states 
(I ■0(0)) = \na) ® \n^) the evolution exibits a cross-phase 
modulation of the form 



\'^{t)) ^ eM~iKnamt)\i,{0)). 



(2) 



We propose an intuitive and straightforward method to 
derive the form of the cross-Kerr interaction in the four- 
level atomic system in a N-configuration for continuous- 
wave fields. By using time-independent non-degenerate 
perturbation theory [isj we find the eigenstates and 
eigenenergies of the time-independent Schodinger equa- 
tion and show that given certain conditions one of these 
eigenstates gives rise to a cross-Kerr interaction. 

The Hilbert space of the four-level atom plus three 
electromagnetic field mode system forms an infinite di- 
mensional tensor product space that is spanned by the 



|2,na-1,nb ,nc> 



|4,na -1,nb+1,nc-1> 




|l,na,nb,nc> 



|3,na-1,nb+1,nc> 



FIG. 1: The four-level atomic system interacting with three 
electromagnetic field modes, and are the Rabi fre- 
quencies and single photon detunings corresponding to each 
electromagnectic field mode. 



set of vectors 



® Wa) ® \nb) (g> \nc)}. 



(3) 



Here i S [1,4] denotes the atomic state and rix S [0,oo) 
describes the Fock state of the corresponding field mode 
fl^. It is assumed that the four- level atom plus three 
electromagnetic field modes form a closed system. This 
enables us to decompose the entire space into four- 
dimensional subspaces that are invariant under the uni- 
tary evolution. That is, given a system in a particular 
state, the evolution will remain in the subspace corre- 
sponding to absorption and creation of single photons 
due to transitions between the four atomic states p^ . 
Omitting the tensor product symbols, the ordered basis 
is given by 



{\l,na,nb,nc), \2,na 
|4, na ~l,nb + I, Tic 



- l,nb,nc), \3,na 
1)} 



l,nb + l,nc), 



(4) 

Here the basis vectors are labeled by the numbers Ux, 
which give the number of photons in the field when 
the atom is in the state Henceforth, we will further 
abbreviate the notation for the basis by refering only to 
the atomic state 



{|1),|2),|3),|4)}. 



(5) 
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We assume that the electromagnetic fields couple to the 
atomic energy levels by the electric-dipole interaction and 
that the rotating-wave approximation can be applied. 
Working in an interaction picture we find that the Hamil- 
tonian can be written in the bare atomic state basis as 



H = h 



ni/2 

ni/2 52 
nc/2 




(6) 



The form of the Hamiltonian is identical to that used in 
semi-classical calculations: the only difference being the 
definition of the Rabi-frequencies which are given by 



2ga\/n^ 



rif, = 2gb\/nb + 1, 
ric = 2gc^/n^. 



(7) 
(8) 
(9) 



Here, is the number of photons in the field mode x 
when the atom is in the state |1) and gx is the coupling 
strength which is proportional to the electric-dipole ma- 
trix elements for each transition. The multi-photon de- 
tunings are given by 



51 = Aa, 

52 = -Ab + Aa, 

53 ^ Ac-Afc + Afl 



(10) 

(11) 
(12) 



To simplify the application of perturbation theory we in- 
troduce a symmetry into lambda subsystem: we assume 
that the fields fia and 0.1, are Raman-resonant with the 
two-photon transition from the state |1) to |3) (^2 — 0). 
This greatly simplifies the determination of the lambda 
subsystem dressed states. In our Raman-resonant case 
the multi-photon detunings reduce to 



Si 



0, 



(13) 



II. 



NONDEGENERATE TIME-INDEPENDENT 
PERTURBATION THEORY 



For the purpose of applying perturbation theory we 
split the system into two parts: the unperturbed system 
described by the Hamiltonian Hq and a weak perturba- 
tion described by V. The unperturbed system consists of 
the four atomic levels coupled by the two fields and 
fif,. This constitutes a three- level lambda-scheme plus a 
fourth uncoupled level |4). 

The weak perturbation V arises due to the third field 
ilc which couples the atomic levels |3) and |4). We as- 
sume that the perturbation strength is dictated by a sin- 
gle multiplicitive parameter e (e = il'^/2). Therefore, the 
Hamiltonian can be written as 



H = Ho + eV. 



(14) 



For the total Hamiltonian given above ([6]) Hq and eV are 



( f7*/2 

Oal2 Si nb/2 

Vll/2 




(15) 



eV 
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e 



(16) 



Before using perturbation theory, it is necessary to know 
the eigenstates of the unperturbed Hamiltonian Hq. As 
mentioned previously, the eigenstates of the lambda- 
system have a particularly simple form for the special 
case of a Raman-resonance between the lower levels [l^ . 
We choose the following normalised eigenbasis: 



(0)n 



1 



^ V al / 



f^a|3)), 

+ ni\3) 
+ ni\3) 

= |4). 

The normalisation constants are given by 

g2 = \na\^ + \nb\^ 

Nl = 2 (g^ + Sl±SiJ5l + G^ 



2A-I2)), 
2A+I2)), 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



and the corresponding eigenenergies (i?o|'/'' 
ij(°)|4°)) = ;,Ai''Vi°^»are 



(0)\ 



AO) 



A 



(0) 
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= A_ 
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= A+ 













G2 



G2 



(23) 
(24) 

(25) 
(26) 



Consider now the dynamical stability of these eigen- 
states. Assuming that the system is weakly coupled to 
the environment, spontaneous emission will relax the sys- 
tem into its most radiatively stable eigenstate, but can 
then be neglected. From inspection we can see that only 

the state \4>f''') is radiatively stable. Moreover, this state 
also exibits no material polarisations between atomic lev- 
els associated with allowed electric-dipole transitions and 
is therefore non-interacting or dark to the applied fields. 
Therefore, the system will relax from any initial mixed 
state into the dark .state \(t)^i'^), or its perturbed counter- 
part. 



We now wish to calculate the approximate eigenener- 
gies and eigenstates of the total Hamiltonian. To find 
approximate solutions we assume that these can be ex- 
panded in a power series of the interaction strength e: 



\K) = Y.^W)- (27) 



i=0 



i=0 



We can expand each term l^n'') of the power series 
in terms of these unperturbed Hamiltonian basis states 
u(o)\. 



s=l 



(28) 



Substituting the expansions ([TT)) into the energy eigen- 
value equation H\<j)n) = E„\(l)n) one finds 

oo oo oo 



i=l 



j=0 k=0 



We assume that this equation is true for terms involv- 
ing each power of e independently. Using (|28p and the 
normalisation condition {4>n\4>'n) = 1; one finds a set of 
coupled equations for the En^ and a^'-*-' terms. These 
equations determine the eigenvalues and eigenvectors of 
the system p^ . 

Let us first consider the perturbation caused to the 
eigenenergies. In the eigenbasis of the unperturbed sys- 
tem the interaction V only has zero elements along the 
main diagonal. Therefore we find that all the eigenen- 
ergies are unchanged to first order. The first non-zero 
contributions are second order corrections given by 



= E 



\{4>f^\V\<t>'n')\ 



(0)\ 



= 1 



E. 



(0) 



This results in the eigenenergy corrections [En 



X 



(2) 



<53(|f^aP + |r!fcP)' 

V2(2A+ + 2(53)' 

iiSi- 63)63 + \na\^ 



63{G^+463{6i~63))- 



(30) 



hx^:^] 



(31) 
(32) 

(33) 

(34) 



For example, the first energy eigenvalue is given to second 
order by 



E^ « hX<f^ + 



'-hX^^^=- 



n\na\^\n. 



A63i\na\^ + W)- 



(35) 



It is possible explain the form of this eigenvalue. Namely, 
the energy correction arises due to the ac-Stark shift of 



101°' 



63 



FIG. 2: The two-level subsystem consisting of the composite 
field f2 = — f2afic/G coupling the ground state [(^i"') to the 
excited state 



the level |3) by the field flc- Assuming that the system 
adiabatically follows this perturbation, then no atomic 
transitions between the unperturbed states will be in- 
duced. The probability of the atom occupying the level 
|3) is If^Qp/G^ if the system is in the state \4>f^). Multi- 
plying this by the ac-Stark shift produced by a weak field 
on that state, — |r2cp/4(53, results in the energy shift 



G2 



4(53 



(36) 



Using ([2T|) . this simple calculation again yields the eigen- 
value previously derived by the second-order perturba- 
tion theory (|35p . 

It should be noted that the eigenvalue correction |3T|) 
fails to provide a good approximation when the detuning 

63 tends towards zero. However, in the summation (j30p 

(2) 

the only term contributing to the energy E^ ' of the dark 
state 14)^1 arises due to the coupling of this state to the 
excited state 104"''). Therefore, to second order in the 

perturbing field Qc, the dark-state It/)'"') behaves as if it 
was the lower level in a two-level subsystem (Fig. [2]). 

We may introduce the ordered basis for the pseudo- 
two-level subsystem 



14"^) = |2) = 



(37) 



The matrix elements of this Hamiltonian are given by 
[HTLs]mn = {m\Ho + eV\n). This gives 



TLS 



= h 



n*/2 



(38) 



n/2 63 

where we define the composite field il — —ila^c/G. For 



4 



a two-level system the lower eigenenergy is given by 



E 




(39) 



This is valid up to second order in fJc, but unlike the non- 
degenerate perturbation theory result this correction also 
holds for all values of detuning 5^ . In the limit of strong 
detuning we again recover the standard perturbative re- 
sult dervied previously ([35|) . 



where we have defined the strength of the cross-Kerr non- 
linearity as 



K ^- 



\9a?\9c 



(46) 



Importantly, the field fif, only modulates the strength of 
the cross-Kerr effect, although the interaction stregth is 
limited by the requirement that \^b\ ^ \^a\, \^c\- 



III. CROSS-KERR INTERACTION 



IV. LIMITATIONS 



We now wish to show that when the system has relaxed 
into the perturbed dark state, a cross-Kerr interaction 
will arise between the fields JIq and Qc- The perturbed 
dark-state is given to first order by 



101 



1 

G 



nb\i)-na\3) 



2^ 



|4) 



(40) 



Essentially, by calculating the approximate eigenenergy 
of the perturbed dark state \4>i) we have determined the 
evolution of the system when relaxed into this eigenstate. 
This evolution is given approximately by 



10(0) «exp(-zi?f)t/ft)|^(0)). 



(41) 



where the intial state is 10(0)) = |0i)- 

To shown that this evolution gives rise to a cross- 
Kerr interaction we place one further limitation upon 
the system. We assume, that the electromagnetic field 
fif, is much stronger than the fields fia and Sic {\^b\ ^ 
I rial, \^c\)- On applying this condition the eigenenergy 
and eigenstate simplify to 



El 



n\na\^\nc\' 

|1) <E) \na) (gi \nb) (g) \nc 



(42) 
(43) 



Recalling again the definition of the Rabi frequencies and 
subsituting these into the energy eigenvalue, we find 



El 



(44) 



Since the eigenstate of the system is now approximately 
equal to the bare state |l,na, nb,nc) we can replace the 
numbers tIq, Ub and in the eigenvalue Ei with the 
corresponding number operators fia, fib and fic. This is 
possible since the state vector \l,na,nb,nc) is an eigen- 
state of the operators h^. with eignevalues rix . 

The evolution of the system is therefore described by 



|7/;(i)) =exp(-iifn„ncO 10^(0)), 



(45) 



It has been assumed throughout that a Raman- 
resonance exists between the lower atomic levels. This 
condition is necessary for the generation of a pure cross- 
Kerr nonlinearity. Moreover, without this constraint the 
system would exibit an additional phase delay and self- 
phase modulation of the field Qa- Due to the narrow 
bandwidth of the EIT window, we expect the system to 
be very sensitive to small changes in two-photon detun- 
ing. For instance, introducing a time-dependence to the 
field fla would invalidate the Raman-resonance condition 
due to the spectral width of the pulse. 

However, since the field Sic need only be strongly de- 
tuned, increasing the spectral width should not have a 
detrimental effect. Each spectral component would expe- 
rience an independent cross-Kerr interaction. Nonethe- 
less, only the field Sla experiences slow-light conditions 
so the interaction time of probe pulses is constrained by 
the group-velocity mismatch. Extending this analysis to 
the non-Raman-resonant regime will be the subject of a 
further publication. 



V. CONCLUSION 

By applying non-degenerate perturbation theory to the 
time-independent Schrodinger equation we have shown 
that the cross-Kerr interaction can arise in an atomic 
four-level system in the N-configuration. In paticular, 
we perturb a Raman-resonant lambda system by intro- 
ducing a fourth atomic level weakly coupled to one of 
the lambda-system ground states. The system relaxes 
into the perturbed dark-state of the lambda subsystem 
and gives rise to the cross-Kerr nonlinearity in the strong 
pump (fib (|Slb| ^ |Sla|, \^c\) and strong detuning limits 
{63 3> Sic). Importantly, this interaction has been shown 
theoretically [l0| and experimentally [llj to be several 
orders of magnitude larger than that produced by other 
contemporary methods. 
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